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ABSTRACT 

We present V, V-I color-magnitude diagrams (CMDs) for three old star clusters in the Large Magellanic Cloud 
(LMC): NGC 1466, NGC 2257 and Hodge 11. Our data extend ~ 3 magnitudes below the main-sequence turnoff, 
allowing us to determine accurate relative ages and the blue straggler frequencies. Based on a differential com- 
parison of the CMDs, any age difference between the three LMC clusters is less than 1.5 Gyr Comparing their 
CMDs to those of M 92 and M 3, the LMC clusters, unless their published metallicities are significantly in error, 
are the same age as the old Galactic globulars. The similar ages to Galactic globulars are shown to be consistent 
with hierarchial clustering models of galaxy formation. The blue straggler frequencies are also similar to those of 
Galactic globular clusters. We derive a true distance modulus to the LMC of (m-M)o = 18.46 ± 0.09 (assuming 
(m-M)o = 14.61 for M 92) using these three LMC clusters. 

Subject headings: blue stragglers — globular clusters: individual (NGC 1466, NGC 2257, Hodge 11) — 
Magellanic Clouds 



1. INTRODUCTION 

Only in the Magellanic Clouds and the Milky Way Galaxy 
' can we directly measure the ages of old star clusters and 
, thereby accurately model the history of cluster formation. In 
■ the Galaxy, after decades of research into dating techniques and 
' collecting data for halo globular clusters (including, very re- 
cently, globulars in the far reaches of the halo using the Hubble 
Space Telescope and WFPC2 [e.g., Harris et al. 1997; Stetson 
et al. 1999)), we are beginning to understand the details of the 
. cluster age distribution. As this story unfolds, a next step is to 
' compare the early formation history of the Galaxy with that in- 
ferred from the same measurement techniques for star clusters 
in the Magellanic Clouds. 

The most important questions to ask are (1) whether the old- 
est Large Magellanic Cloud (LMC) clusters are as old as the 
oldest Galactic globular clusters (GGCs), and (2) whether there 
is an age spread among the most metal-poor LMC clusters. 
The answer to the first question tells us about the epoch of ini- 
tial cluster formation for galaxies of different mass and Hub- 
ble type. The age spread among the old clusters is one tracer 
of the early star formation history of the galaxy, and indicates 
timescales important for collapsing or merging gas clouds. Al- 
though the complete age distribution of the GGCs is not yet 
known, the evidence to date suggests that the majority of the 
Milky Way's halo clusters formed at the same time, with a small 
fraction of clusters that formed up to ^ 4 Gyrs later (Stetson, 
VandenBerg & Bolte 1996; although see Sarajedini, Chaboyer 
& Demarque 1997 for a different view). What fraction of the 
co-eval clusters formed in a single parent body or in smaller 
structures that later merged in not clear. The present results for 
the LMC also suggest a mostly co-eval population of old clus- 
ters (Brocato et al. 1996; Olsen et al. 1998). 

' Guest User, Canadian Astronomy Data Centre, which is operated by the H( 
'Bolte, P.I. Proposal Number 05897 



The first studies to attempt to measure ages for the old 
clusters in the LMC appeared in the mid 1980s. Studies of 
Hodge 11 (Stryker et al. 1984; Andersen, Blecha, & Walker 
1984), NGC 2257 (Stryker 1983; Hesser et al. 1984), and sev- 
eral clusters in the series of articles by Walker (referenced in 
Walker, 1993) led the way in studies of these sorts. However the 
combination of large distance to the LMC and problems with 
field-star contamination made even relative age determinations 
with precision of a few Gyr or better difficult. Because of the 
superb resolution of the Hubble Space Telescope (HST), color- 
magnitude diagrams (CMDs) can now be measured in the clus- 
ter cores, which substantially reduces the impact of field star 
contamination. With a modest effort using HST, stars three to 
four magnitudes below the main-sequence turnoff can be mea- 
sured with good accuracy. 

We have WFPC2 images in F555W and F814W of seven 
old LMC clusters selected to have RR Lyrae stars or very blue 
horizontal-branch morphology (see Table 1). In this paper we 
present color-magnitude diagrams for NGC 1466, NGC 2257 
and Hodge 11, estimate their ages relative to GGCs, and deter- 
mine their blue straggler specific frequencies. 

2. OBSERVATIONS AND DATA REDUCTION 

All of our clusters were observed in Cycle 5.' Each cluster 
was imaged for 2 x 260s and 3 x 1000s in F555W John- 
son V) and 2 x 260s and 4 x 1000s in F814W (- Cousins /). 
The exposures were dithered to reduce the effect of undersam- 
pling which is aggravated by subpixel quantum-efficiency vari- 
ations. For NGC 1466, the cluster was centered on the PC. Both 
NGC 2257 and Hodge 1 1 were centered on WF3. Figures la - 
c show the mosaicked images of our three clusters. 

For our final reduction, we obtained recaUbrated images from 

■g Institute of Astrophysics, National Research Council of Canada. 
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Table 1 
LMC Cluster Data 



Cluster 


[Fe/H] 


Mv 


Proj. Radius 


RR Lyrae? 


Hodge 1 1 


-2.06 


-6.99 


4.7° 


no 


NGC 1466 


-1.85 


-7.89 


8.4° 


yes 


NGC 1841 


-2.11 


-7.86 


8.4° 


yes 


NGC 1786 


-1.87 


-7.88 


2.5° 


7 


NGC 2210 


-1.97 


-8.09 


4.4° 


yes 


NGC 2257 


-1.80 


-6.91 


8.4° 


yes 


Reticulum 


-1.71 


-5.96 


11.4° 


yes 



Note. — Data from Suntzeff et al. (1992) 



Fig. 1 



— Mosaicked V-band images of the NGC 1466 WFPC2 field (top-left), the NGC 2257 WFPC2 field (top-right) and the Hodge 1 1 WFPC2 field (bottom-left) 
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the Canadian Astronomy Data Centre that had have been re- 
reduced using the most appropriate biases, darks, flats and bad- 
pixel masks. These calibration frames were obtained closer in 
time to the object frames than the ones originally used in the 
standard HST pipeline. We then masked the vignetted portions 
of the chips and corrected the pixel values for changes in effec- 
tive pixel area over the chips. 

For flie photometry, we used the ALLSTAR/ALLFRAME 
suite of programs which do point-spread-function (PSF) fitting 
to determine magnitudes (Stetson 1987, 1992, 1994). The PSF 
varies spatially across each chip and temporally as the HST fo- 
cus drifts. We created a PSF by combining Tiny Tim (Krist 
1995) PSFs and PSFs generated by DAOPHOT from our im- 
ages. A DAOPHOT PSF for the inner three pixels of the model 
PSF was created for each chip using the images that had the 
best distribution of stars on that chip. Because we lacked iso- 
lated, bright stars to determine the wings of the PSF, we used 
the Tiny Tim PSF beyond 3 pixels. Judging by the widths of 
the horizontal and red giant branches and by the lack of trends 
of aperture correction with position, our PSFs did a reasonable, 
although not perfect, job of mapping the spatial variations in the 
PSF. We will return to this when discussing our fiducial lines. 

We set the weighting parameters "profile error" and "percent 
error" to zero in ALLSTAR and ALLFRAME, as suggested by 
Cool & King (1995). This eliminates the deviation of the pixels 
from the PSF and errors in the flatfielding from the weighting 
scheme. ALLSTAR and ALLFRAME essentially did intensity- 
weighted aperture photometry of the inner 1.75-2 pixels, then 
added a correction for the amount of light in the wings appro- 
priate for the PSF at the star's position. 

In detail, for each cluster we ran FIND and ALLSTAR on 
each frame separately. Next, we matched stars between frames 
using DAOMASTER (Stetson 1993) only including stars on 
the master list which appeared on at least three frames. This 
successfully removed cosmic rays from the Ust. This master 
Ust was input into ALLFRAME which uses one star list for all 
frames, but determines magnitudes for each frame individually. 
Thus we have five measured V magnitudes and six I magnitudes 
for most of the objects in each cluster. 

The ALLFRAME magnitudes required several corrections 
before we could calibrate the data using the transformation 
equations of Holtzman et al. (1995). We first needed to correct 
for the difference between our PSF magnitudes and apertures 
of 0".5, which is the aperture size used by Holtzman et al. in 
deriving their transformation equations. To do this, we found 
cosmic ray hits on our images by comparing each individual 
image with a cosmic -ray-cleaned image created by taking the 
median of a registered stack of frames. Pixels in each frame 
which deviated by more than 2cr from the value in the medi- 
aned image were flagged and the IRAF task FIXPIX was used 
to interpolate over them. We next selected a set of fairly iso- 
lated stars on each frame and subtracted all other stars. Stars 
that had more than two flagged pixels within 0''5 or any pix- 
els flagged within three pixels (~ two PSF fitting radii) of their 
centers were rejected. For each frame, the mean difference be- 
tween the PSF and 0''5 aperture magnitudes for the culled list 
of stars was determined. This value was the aperture correction 
and was added to all the PSF magnitudes for that frame. 

For the short exposures on the PC chip for NGC 2257 and 
Hodge 11, the difference between aperture and PSF-based mag- 
nitudes was magnitude-dependent. This is due to the WFPC2 
charge-transfer efficiency (CTE) problems, since there were not 
enough photons in the wings of the faint stars to fiU up the 



charge traps. For these two clusters the PC has only a small 
fraction of the total cluster stars, so we did not use it. For 
NGC 1466, the PC had a much higher background and the aper- 
ture corrections were constant with magnitude. We also note 
that we used sky annuli of 2"0 to 2"5 , instead of the 4" and 
6" that Holtzman et al. used. This reduced the effects of badly 
subtracted neighbors in crowded regions and uneven sky bright- 
ness across the cluster face. However, there is a difference be- 
tween using the closer sky apertures and the more distant ones 
due to contributions from the wings of the stars. Based on our 
model PSFs, this difference is only 0.001 mags - much smaller 
than our errors in the aperture corrections (see below). 

Within 0"5, most of stars were affected by residuals from 
subtracted neighbors and cosmic rays. The RMS scatter in the 
aperture corrections (aperture - PSF magnitudes) for an indi- 
vidual frame ranged from ~ 0.02 to ~ 0.04 mag. To see what 
random errors this observational scatter introduced into our ze- 
ropoints, we used DAOMASTER to calculate the average mag- 
nitude offset for each frame among all the stars common to all 
frames in each filter + field combination. After accounting for 
the expected offset due to different exposure times, any remain- 
ing offset is a result of errors in our aperture correction. We are 
assuming here that the PSF of WFPC2 did not change in the 
amount of time it took to take images of one cluster (~ 4 hours). 
We fit a Gaussian to the distribution of all the offsets and found 
our 2cr errors to be 0.014 mags. Since DAOMASTER removes 
these offsets when averaging magnitudes, essentially bringing 
all images to the zeropoint of one frame, this is the error for 
our averaged magnitude zeropoints as well. We note that this 
error does not include any systematic error that may result if 
0"5 apertures do not include a PSF-independent fraction of the 
light. 

After adding the aperture corrections, we corrected for the 
CTE problem in the y-direction, using the prescription of Whit- 
more & Heyer (1997) for 0"5 apertures. Our PSFs are not 
good enough to allow accurate interpolation of saturated stars, 
so these stars were eliminated. We averaged our five F555W 
measurements and six F814W measurements using DAOMAS- 
TER and kept only the stars with errors < 0.08 mag in both 
filters. Finally, we caUbrated the data using the Holtzman et al. 
equations. We list the photometry and positions for all our stars 
in Tables 2a-c. We used the first V exposure to determine the 
Usted X and Y coordinates. (Tables 2a-c available from the first 
author) 

In addition to the random errors caused by scatter in the 
aperture corrections (0.014 mag), there is added uncertainty 
from the Holtzman et al. zero points (~ 0.03 mag), and from 
the Whitmore & Heyer corrections (~ 0.02 mag) (see Stetson 
1998). Adding these errors in quadrature leads to an overaU 
error of 0.04, most of it systematic, in our zeropoints. The sys- 
tematic error in V-Zis more likely to to be ^ 0.02, because the 
uncertainties in the Whitmore & Heyer corrections and proba- 
bly in the Holtzman zero points are correlated between V and 
I. 

3. COMPARISON WITH PREVIOUS PHOTOMETRY 

Each of these clusters has at least one CCD-based photomet- 
ric study that we can compare with ours. Our general procedure 
for making comparisons was to find stars in common between 
our study and the previous ones and eliminate all stars from the 
joint sample where we detected a companion within 0"5 in our 
HST images. We also exclude stars within 1 magnitude of the 
detection limit of the ground-based sample. In each case our 
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Fig. 2. — Comparison of Walker's and our photometry for NGC 1466. A line is drawn at Amag = for reference. 



Table 3 
Photometry Comparisons 



Cluster 


Author 


Chip 


Filter 




# Stars 


NGC 1466 


Walker 1992b 


WF2 


V 


-0.007 ±0.018 


32 


NGC 1466 


Walker 1992b 


WF3 


V 


0.009 ±0.020 


17 


NGC 1466 


Walker 1992b 


WF4 


V 


-0.008 ±0.016 


16 


Hodge 1 1 


Walker 1993 


WF2 


V 


0.016±0.011 


80 


Hodge 11 


Walker 1993 


WF4 


V 


0.039 ±0.014 


53 


Hodge 11 


Walker 1993 


WF2 


I 


-0.007 ±0.026 


87 


Hodge 11 


Walker 1993 


WF4 


I 


-0.007 ±0.021 


54 


Hodge 11 


MighellefaZ. 1996 


WF3 


V 


0.051 ±0.005 


144 


Hodge 1 1 


Mighellefa/. 1996 


WF4 


V 


0.032 ±0.007 


46 


NGC 2257 


Walker 1989 


WF2 


V 


-0.017 ±0.008 


26 


NGC 2257 


Walker 1989 


WF3 


V 


-0.027 ±0.013 


83 


NGC 2257 


Walker 1989 


WF4 


V 


-0.002 ±0.008 


16 


NGC 2257 


Testa effl/. 1995 


WF2 


V 


-0.047 ±0.008 


102 


NGC 2257 


Testae/ a/. 1995 


WF3 


V 


-0.012 ±0.005 


100 


NGC 2257 


Testae? a/. 1995 


WF4 


V 


-0.046 ±0.008 


109 



calculated offsets and errors are tabulated in Table 3. stars with V < 23.2 mag and further than 13 "6 from the cluster 

For NGC 1466, Walker (1992b) pubUshed B, V values for center. Figure 2 shows the differences between photometries 
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Comparison of Walker's and our photometry for Hodge 11 . (a) V magnitudes (b) I magnitudes. A line is drawn at Amag = for reference. 
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Fig. 4. — Comparison of Mighell etal.'s (1996) and our photometry for Hodge 11. We plot the M96 stars fovmd on our (a) WF3 and (b) WF4. A line is drawn at 
Amag = for reference. 
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Fig. 5. — Comparison of our photometry and (a) Walker's photometry (b) Testa et al.'s photometry for NGC 2257. A line is drawn at Amag = for reference. 



for the three WF chips. (The PC had very few stars with mea- 
surements by Walker.) The stars in common between the two 
studies are mostly at the faint end of Walker's study where he 
did not quote errors. So instead of finding a weighted mean, we 
found the median Amag. Our estimate of a comes from half of 
the difference between the value > 16% of the Amags and the 
value > 84% of the Amags. This estimate of a is then reduced 
by square root of the number of measurements, to estimate the 
standard error of the mean. 

For Hodge 11, Walker (1993) imaged in V and I. Because of 
the crowding from the ground. Walker excluded all stars closer 
to the cluster center than 40", so no stars from WF3 were in- 
cluded in the comparison. The V and / magnitude comparisons 
(Figures 3a-b) show a large scatter, but no zeropoint offset or 
trends with magnitude. Mighell et al. (1996) published B, 
V data on Hodge 1 1 obtained as part of a snapshot survey of 
Magellanic Cloud clusters. In Figures 4a-b, we plot the differ- 
ence between our WF3 and WF4 data and matching stars in the 



Mighell et al. photometry. We calculated a weighted average 
of difference, using 3a clipping. 

NGC 2257 is the sparsest cluster of the three and has ground- 
based photometry in the cluster center. Both Walker (1989) and 
Testa et al. (1995) have B, V data. Testa et al. used Walker's 
photometry to calibrate their data, but the Testa et al. data ex- 
tend two magnitudes fainter for stars further than 40" from the 
cluster center Qualitatively, Figures 5a-b show that our pho- 
tometry is in reasonable agreement with previous efforts. Our 
zeropoint error estimate of 0.04 mag seems reasonable in light 
of the comparison with previous photometry. We note that rel- 
ative age determinations and blue straggler statistics are unaf- 
fected by absolute scale concerns, although the reddenings and 
distances determined will vary for different calibrations. We 
see no signs of non-linearity or color transformation problems 
that could affect the results discussed in this paper. 

4. COLOR-MAGNrrUDE DIAGRAMS 
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Fig. 6. — NGC 1466 color-magnitude diagrams for the individual chips of WFPC2. The cluster is centered on the PCI. 
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Fig. 7. — NGC 2257 color-magnitude diagrams for the individual chips of WFPC2. The cluster is centered on the WF3. 



In Figures 6-8, we present the calibrated CMDs for the clus- 
ters. Zero-point differences cause slight offsets between the 
CMDs derived from the various chips within a given cluster. 
The largest shifts are ~ 0.025 mag in color, which can almost 
all be ascribed to the 0.015 mag random errors in the aperture 
corrections in V and /, with some error due to the Holtzman 
gain ratios and Whitmore & Heyer (1997) CTE corrections. 
In all cases we have well-defined cluster sequences from near 
the tip of the RGB down to ~ 3 mag fainter than the main- 
sequence turnoff (MSTO). We determined the fiducials for the 
stars on each chip using an algorithm described in Sandquist et 
al. (1996). The main-sequence fiducial was determined by bin- 
ning the data in magnitude and finding the mode in color. The 
red-giant-branch (RGB) fiducial was found by using the mean 



color of magnitude bins. The horizontal branch (HB) and sub- 
giant branch (SGB) could not be reliably fit by such methods, 
the HB because of the large color dispersion of the RRLyr and 
the SGB because it was neither vertical or horizontal. In these 
two regions we found points by eye, but the spread in magni- 
tude in these areas is small, and this procedure should not add 
large errors. For chips with a smaller number of stars, the bins 
used for finding the mode and mean were larger, and the fidu- 
cials are noisier. In Figures 9a-c, we plot the fiducials from each 
chip for each cluster For NGC 2257 and Hodge 11, the clus- 
ters which have their centers on WF3, we used a sample with 
r > 20" from the cluster center to determine the WF3 fiducial. 
For NGC 1466, we used stars with r > 9"2 from the cluster 
center for the PCI fiducial. 
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Table 4 
Cluster Fiducial Points 



NGC 1466 NGC 2257 Hodge 1 1 

V V-I V V-I V V-I 
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Fig. 8. — Hodge 11 color-magnitude diagrams for the individual chips of WFPC2. The cluster is centered on the WF3. 



In our plot for NGC 2257, we include the fiducial determined 
from the inner 20" on WF3. It is offset by a constant color from 
the fiducial determined from the stars with r > 20" . We believe 
this is a manifestation of our less-than-perfect PSFs, which do 
not fully account for the changes in the PSF over the face of 
the chip. However, examination of Figures 9a-c shows that this 
is not a major concern, since the fiducials from each chip and 
each region look very similar. 

For our final fiducials, we adopt the fiducial of the chip con- 
taining the cluster center (Table 4). We also calculated the 
shifts needed to compensate for chip-to-chip zero-point offsets 
by matching the SGBs and RGBs of the individual chip fidu- 
cials. These shifts were applied to the data to create a master 
CMD of all stars on all chips to use when finding blue strag- 



glers. 

5. RELATIVE AGES 

5.1. Metallichies 

In order to inter-compare the CMDs of the LMC clusters and 
to compare them with those for the GGCs, we need to know 
their metallicities. We will refer all measurements to the Zinn 
& West (1984; ZW84) scale. The accuracy of this scale has 
been questioned by Carretta & Gratton (1997). They deter- 
mined [Fe/H] for individual stars in 24 GGCs using a homo- 
geneous analysis of equivalent widths of Fe lines from high- 
resolution spectra. Their correlation with the ZW84 metal- 
licities is significantly non-linear, especially for intermediate 
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Fig. 9.— 



The fiducials determined separately from the data for each chip, (a) NGC 1466 (b) NGC 2257 (c) Hodge 11. 



JOHNSON ET AL. 



9 



metallicities around -1.5. Rutledge et al. (1997) compared 
metallicities based on W' , the reduced equivalent width of Ca 
II triplet lines, with ZW84 and Carretta & Gratton. This es- 
timate of the metallicity also was related non-linearly to the 
ZW84 [Fe/H]'s, but linearly to the high-resolution results. This 
result is somewhat ambiguous because of the dependence ofW' 
on log g and [Ca/Fe], but shows again that the relative, in ad- 
dition to the absolute, scale of GGC metallicities has not been 
settled. However, we are interested in a small enough range 
that the ZW84 scale can be considered linear Also, for relative 
age studies, sensitivity to [Fe/H] uncertainties decreases with 
decreasing Z and our results are not sensitive to [Fe/H] errors 
< 0.3 dex. 

There are only a few studies of the abundances for the old 
LMC clusters. Cowley & Hartwick (1982) used low-resolution 
spectra to measure the G-band strength, Ca II break and the 
average strength of groups of weak Fe lines. They calibrated 
indices against several GGCs and used the Zinn (1980a) scale, 
which is almost identical to the ZW84 scale. Olszewski et 
al. (1991) measured the strength of the Ca IR Triplet lines 
for a large number of LMC clusters, including NGC 1466 and 
Hodge 11. They derive [Fe/H] via comparison with the Ca 
strength in stars in the GGCs NGC 288, NGC 1851, M 79 and 
NGC 7099 whose [Fe/H] were from ZW84. 

For three stars in NGC 1466, Cowley & Hartwick found 
[Fe/H]=-2.0 ± 0.2. Olszewski a/. (1991) derived abun- 
dances for two stars, LE4 and LE3, in NGC 1466 and found 
that LE4 had a substantially lower [Fe/H] (-2.48) than LE3 
(-1.85). Walker (1992b) argued that since the apparent magni- 
tude of LE4 placed it significantly above the cluster sequences, 
it probably had a close, bluer companion which weakened the 
Ca II equivalent widths. Checking our HST images, we find 
LE3 is reasonably well isolated, while LE4 is indeed blended 
with at least three other objects that are closer than 0"5. Walker 
also estimated [Fe/H] via several methods involving the proper- 
ties of NGC 1466 RR-Lyrae stars and the de-reddened colors of 
various [Fe/H]-sensitive points in the cluster CMD. He found 
a consistent value of -1.82 ± 0.04 from all of his methods if 
E(fi-y ) in the direction of NGC 1466 is ~ 0.08mag. 

Cowley & Hartwick measured three stars in Hodge 1 1 and 
derived a mean value of [Fe/H]= -2.1 ±0.2, consistent with 
Olszewski et al.'s average from two stars of -2.06 ±0.2. In 
this case, both stars had no companions within 0"5 on our HST 
frames. Based on the de-reddened color of the giant branch in 
V-I , Walker (1993), using the calibration of Da Costa & Ar- 
mandroff (1990), derives [Fe/H]= -2.0 ± 0.2 for E(V-I )= 0. 1 . 

For NGC 2257, Cowley & Hartwick measured a mean 
[Fe/H]= -1.8 ±0.3 from five stars. Testa et al. measured the 
metalUcity indicator {B-V)o,g (Sandage & Smith 1966) from 
their CMD which leads to a value of -1 .86 on the ZW84 scale. 
Walker (1989) used the period-A^ diagram of the RR Lyraes to 
estimate the metalhcity as -1.8 ±0.1. The agreement between 
studies for each cluster is excellent. 

We will adopt [Fe/H]=-2.05 for Hodge 11 and -1.85 for 
NGC 1466 and NGC 2257. These measurements are on the 
ZW84 scale, but in this |Fe/H| regime there is little disagree- 
ment between ZW84 and any of the more recent scales. The 
GGCs that we use for comparison are M 92 (-2.25) and M 3 
(-1.66) (ZW84). 

It is of potential concern that the LMC clusters could have 
different abundance ratios, such as [O/Fe], [C/Fe] and [a/Fe], 
than their Galactic counterparts. However, Cowley & Hartwick 
noted no difference in the metallicities derived by comparing 



the strength of weak Fe lines to those of Galactic clusters and 
the metallicities found by looking at the Ca II break and the G- 
band, which are dominated by the abundance of the a elements. 
In the absence of any evidence to the contrary, we will assume 
the same trends in abundance ratios as are in GGCs. 

5.2. Ages: CMD comparisons 

Figures 10-12 show the fiducials for the LMC clusters com- 
pared to V, y-7 fiducials for M 92 and M 3 (Johnson & Bolte 

1998). 

We accounted for the differences in distance and small dif- 
ferences in abundance by registering the HE magnitudes at the 
blue edge of the RR Lyr gap so that AVhb = -Q.2QA[Fe/H]. 
We accounted for the differences in reddening by making the 
relative colors of the RGB agree with Bergbusch & Vanden- 
Berg (1992) theoretical isochrones of the correct metalhcity. 
In Figure 10a, the NGC 1466 fiducial line fits between M 92 
and M 3 from the RGB (of course this is guaranteed by our 
procedure for correcting for E(V-I) on the upper RGB) down 
through the MS as expected for a cluster with a metallicity inter- 
mediate between M 92 and M 3 and the same age as these two 
clusters. In Figure 10b, we show how well we can determine 
an age difference from such plots. We have taken Bergbusch 
& VandenBerg (1992) (BV92) isochrones for [Fe/H]=-2.26, 
-1.78 and -1.66. On this plot, we represent M 92 and M 3 
by 12 Gyr isochrones of [Fe/H]=-2.26 and -1.66. Then we 
plot [Fe/H]=-1.78 isochrones for 13, 12, 11, 10 and 9 Gyrs. 
The differences in the turnoff region are striking and show that 
NGC 1466 has the same age as M 3 and M 92 to < 1.5 Gyr. 
For NGC 2257 and Hodge 11, Figures 1 1-12 show the compar- 
isons, and we reach similar conclusions. 

5.3. Ages: Color-Difference Method 

We can quantify the age differences between the GGC and 
the LMC clusters and among the LMC clusters themselves by 
measuring the color difference between the MSTO and the base 
of the RGB (VandenBerg, Bolte, & Stetson 1990 (VBS); Sara- 
jedini & Demarque 1990). In practice, we do this by shifting 
the cluster fiducials horizontally until the colors of the MSTO 
match and vertically until a magnitude reference point matches. 
Since the MSTO does not have a well-defined magnitude, VBS 
suggested using the point on the MS that is 0.05 redder than the 
MSTO color (V+o.os). All other things held constant, an older 
cluster will have a shorter SGB and, after MSTO registration, 
its RGB will lie to the blue of that for a younger cluster. Unfor- 
tunately, the length of the SGB also changes with metallicity. 
While this is a relatively small effect in B-V (VBS), it is larger 
when V-I colors are used. In V-I , the [a/H] is the crucial abun- 
dance to know (VandenBerg, private communication) and in all 
cases we have a good relative metallicities from Ca lines. To 
determine the magnitude of the metallicity effect requires the 
use of theoretical isochrones, and in particular, theoretical col- 
ors. Since colors are among the most uncertain quantities pre- 
dicted by theory, we will use two separate sets of isochrones 
to evaluate our results and their dependence on the choice of 
isochrones. 

Our LMC fiducials are in many ways an ideal dataset to use 
with this method. They were all taken with the same instru- 
ment and have the same data quality. They have low metallici- 
ties within 0.2 dex of each other, so the effects on the colors are 
minimized. We measured the MSTO color by fitting a parabola 
to the individual points in that region. We then interpolated 
between our fiducial points on the MS to determine Vb.os- In 
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V-I V-I^ 

Fig. 10. — (a) NGC 1466's fiducial compared with M 92 and M 3. We have two estimates of the error in the NGC 1466 points. First, we determined the average 
deviation of the PC fiducial from the WF fiducials in Figure 9a. Second, we estimate the random error in the color of each bin to be 0.005 magnitudes. We adopted 
the larger number as our error. The V errors in the HB magnitude are smaller than the boxes used. However, the M 92 and M 3 HBs have larger, but unknown, errors 
because of the relatively few stars that define them (see Johnson & Bolte 1998 for more discussion). (In particular, the blue end of M3's HB is defined by one star.) 
(b) Bergbusch & Vandenberg (1992) isochrones compared in a similar manner. Isochrones of 12 Gyrs are plotted for [Fe/H]=-2.26, -1.78 and -1.66. To illustrate the 
effect of age, [Fe/H]=-1.78 isochrones are also plotted for 9, 10, 1 1, and 13 Gyrs. ([Fe/H]=-1.78 is more metal-rich than the values (c.f. Table 1) for our clusters, so 
the correct isochrones would be bluer and brighter.) 




v-i 



Fig. 11. — NGC 2257's fiducial compared with M 92 and M 3. See the caption for Figure 10 for a discussion of the errors. 




v-i 



Fig. 12. — Hodge ll's fiducial compared with M 92 and M 3. See the caption for Figure 10 for a discussion of the errors. 

Figure 13 we have registered the LMC cluster fiducials using these two points. The sequences for NGC 1466 and NGC 2257 
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Fig. 13. — NGC 1466, NGC 2257 and Hodge 1 1 fiducials compared using the color-difference method. The remarkable agreement in the position of their RGBs, 
especially for NGC 1466 and NGC 2257, argues against a large age difference between the clusters. 



are indistinguishable. Hodge 11 's RGB lies to the blue of the 
other two. With the color-difference method in V-I , the more 
metal-poor Hodge 11 should be redder than NGC 1466 and 
NGC 2257, so this indicates that Hodge 1 1 is somewhat older 
than the other two. As shown below, the derived age difference 
is not large. 

Next, we quantified the age difference between the LMC 
clusters and the GGCs based on any shifts of the postions of 
the registered RGBs. We first fit the GGCs RGBs from the base 
to the magnitude of the HB with a straight line. Next, we cal- 
culated the weighted average of the color offset between these 
lines and the individual RGB stars below the HB of each of our 
LMC clusters. In Figure 14, we show an example of the method 
where we have registered M 92 and NGC 2257 using the two 
fiducial points. Based on VBS, the mean offset for a cluster 
gives the age difference (compared to M 92) and the standard 
error of the mean offset gives an estimate of our 6 age precision. 
The mean offsets for each cluster are recorded in Tables 5a-b as 
6iV-I). 

While the standard error of the mean represents the error in 
our measured shifts caused by scatter on the RGBs, it does not 
take into account errors in our choice of V-I msto and V+o.o5- 
The error in V+o.os caused by scatter on the MS has only a small 
effect on the color shift of the RGB since the RGB is almost 
vertical. A shift of 0.05 mag in y+0.05 results in a color shift of 
0.0025 mag. The more significant error is in V-I mSTO. This 
color could change by an amount of order V-I ^ 0.003 mag if 
we made plausible alterations in the fitting range and the ac- 
ceptable color error for our parabolic fit. We conservatively 
estimate our observational error in V-I msto as 0.005 mag. Un- 
fortunately, an error in the color also results in an error in our 
choice of V+o.os. These two errors work in the same direction. 
For example, if our MSTO color is too red, our V+o.os will be 
too faint. We will then shift our fiducial too far to the blue to 
get our red MSTO color to match and too bright (and hence too 
blue) to get our y+0.05 to match. Taking into consideration the 
slopes of the MS and the RGB, we find that our 0.005 mag er- 



ror in color results in an additional 0.0035 mag shift in the RGB 
color due to the correlation of these errors. Our total random 
error is therefore 0.009 mag. Finally, we note that since we cal- 
culate the difference for our three LMC clusters relative to one 
GGC at a time, the errors in the GGC points will result in a sys- 
tematic shift. Comparing the results obtained for M 92 in Table 
5a with those for M 3 in Table 5b will provide an estimate of 
the size of that effect. 

To convert our measured S(V-I) and errors into age dif- 
ferences and uncertainties, we used two sets of isochrones: 
Bergbusch & VandenBerg (1992) (BV92) and the new Yale 
Isochrones (Chaboyer et al. 1995) (C95). We determined the 
predicted shifts in RGB color as functions of metallicity and 
age after registering the isochrones in the same manner as our 
data. For -1 .66 > [Fe/H] > -2.26, we found that the color shift 
predicted by BV92 depends on the metaUicities by the follow- 
ing relation: 

5{V-I ) = -0.069(A[Fe///]). (1) 

Next, comparing sets of isochrones with the same metallicity, 
but with ages ranging from 9 to 1 3 Gyrs, we found that the RGB 
shift is related to the age difference (in Gyr) by 

S{V-I ) = -0.021(AA^e). (2) 

We first removed the known shift due to metallicity differ- 
ences from our results in Tables 5a-b. The remaining shift is 
due to age differences. Using our &ge,-5{V-I ) relation above, 
we calculated the age differences relative to the GGCs and in- 
cluded those values in Tables 5a-b. Our age-(5(y-/) relation is 
valid only for ages between 9 and 13 Gyrs. Our use of younger 
isochrones to calibrate the RGB shifts is not crucial for our dis- 
cussion of relative ages, but we note that the choice of an older 
age range would increase the implied age difference for a given 
5{V-I ). However, the errors would also increase, and we would 
be left with the same conclusion. 

We performed the same procedure using the C95 isochrones 
with two exceptions. First, the C95 isochrones are not a- 
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Fig. 14. — The color-difference method applied to NGC 2257 and M 92. We plotted the individual stars for each cluster. The solid trapezoid shows the M 92 stars 
used to calculate the straight line fit to the RGB. The dashed trape zoid encloses the NGC 2257 stars used to calculate the RGB shift relative to the M 92 Une. 



Table 5a 
6iV-I) Quantities 



Cluster 


{V-I)to 


V+fl.05 


6(V-I) 


A Age (Gyrs) 


A Age (Gyrs) 










(BV92) 


(C95) 


M92 


0.553 ±0.008 


19.45 ±0.05 


0.000 ±0.014 


0.0 ±0.7 


0.0 ±0.7 


NGC 1466 


0.623 ±0.005 


23 .73 ±0.05 


-0.012 ±0.009 


-0.8 ±0.4 


-1.0 ±0.5 


NGC 2257 


0.610 ±0.005 


23.38 ±0.05 


-0.008 ±0.009 


-1.0 ±0.4 


-1.2±0.5 


Hodge 11 


0.628 ±0.005 


23.58 ±0.05 


-0.021 ±0.009 


+0.3 ±0.4 


+0.3 ±0.5 



Table 5b 
5iV-I) Quantities 



Cluster 


(V-I)to 


V+0.05 


S(V-I) 


A Age (Gyrs) 
(BV92) 


A Age (Gyrs) 
(C95) 


M3 

NGC 1466 
NGC 2257 
Hodge 11 


0.595 ±0.005 
0.623 ±0.005 
0.610 ±0.005 
0.628 ±0.005 


19.92 ±0.05 

23.73 ±0.05 
23.38 ±0.05 
23.58 ±0.05 


0.000 ±0.009 
0.008 ±0.009 
0.011 ±0.009 
-0.002 ±0.009 


0.00 ±0.4 
+0.24 ±0.4 
+0.10±0.4 
+1.28 ±0.4 


0.00 ±0.5 
+0.36 ±0.5 
+0.21 ±0.5 
+1.56 ±0.5 
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enhanced, so we followed their recommendation and deter- 
mined the appropriate solar-scaled isochrones by using the for- 
mulation of Solaris, Chieffi, & Straniero (1993) and choosing 
the relative a abundances to be 4 x solar. Second, we found 
a noticeable trend in the steepness of the slope of the age- 
(5y-/ relation depending on the metallicity of the isochrones. 
However, this affected our ages only at the level of 0.05 Gyrs, 
so we adopted the slope determined from the [Fe/H]= -1.66 
isochrones. Our equations for 5{V-I) versus metallicity and 
age are 

^(y-/) = -0.080(A[Fe///]), and (3) 

5{V-I) = -OmO{Age) (4) 

The C95 isochrones have a steeper dependence on metallicity, 
but the overall agreement between BV92 and C95 is quite good. 
The cahbration of the color-difference method will need to up- 
dated again as newer isochrones become available, but there 
is no evidence that the theoretical calibration will be an over- 
whelming source of error. Based on the SiV-I) method, we 
conclude that the LMC clusters are the same age as the GGCs 
to within 1.5 Gyrs. 

5.4. Ages: Veto 

The relative ages of GCs can also be calculated by using the 
difference in V magnitude between the HE and the MSTO (the 
Ay method) (Sandage 1982). Because luminosities can be pre- 
dicted more reUably than temperatures, this method has been 
often put forward as the most robust of the relative age esti- 
mators. However, it is relatively more uncertain for clusters 
with very blue HB morphologies such as Hodge 1 1 for at least 
two reasons. First, it is difficult to assign an apparent bright- 
ness to the HB and to correct for HB-star evolution in blue-HB 
clusters. Second, as pointed out recently by Sweigart (1997) 
and Langer, Bolte & Sandquist (1999), there is reason to sus- 
pect that some blue HB stars have had helium mixed into their 
envelopes resulting in higher HB luminosity and, for clusters 
containing He-mixed stars, an artificially inflated AV value. In 
addition, our M 3 and M 92 data do not have many HB stars, 
which limits the precision of our HB magnitudes. Finally, de- 
termining the magnitude of the MSTO is uncertain, since the 
color-magnitude diagram is vertical at that point. 

Chaboyer et al. (1996a) proposed a modification that moved 
the lower reference point to the subgiant branch, 0.05 mag red- 
der than the MSTO color (Veto)- This reduces the errors, since 
the CMD cluster sequence is sloped at that point and a fiducial 
point can be determined with more confidence. To find Vbto> 
we used the mean magnitude of the stars in a box on the sub- 
giant branch centered on our fiducial at (y-/)Msro + 0.05. The 
box was 0.04 mag wide in color and 0.1 5 mag wide in V. The 
standard error of the mean for those points within the box pro- 
vides one estimate of the error, but probably underestimates it 
because there is additional confusion at the blue end of the SGB 
due to the scatter of stars from the MS to the red. This makes the 
SGB wider in magnitude and makes Veto dependent on our ex- 
act box size in vertical direction. Our total errors are therefore 
closer to ~ 0.03 mag. For the HB magnitude of NGC 1466 and 
NGC 2257, we found the average and the standard error of the 
mean for our HB stars at the blue end of the RRLyr gap. This is 
the method by which we found the HB magnitude for our GGC 
data, and avoids the problem of determining < My (/?/?) > with 
our limited time resolution. We used the HB of NGC 2257 



as a template to extend Hodge 11 's extremely blue horizontal 
branch to the red. We ignored the two reddest points of the 
Hodge 1 1 HB as Ukely being either evolved stars or field in- 
terlopers. The Hodge 1 1 HB magnitude derived in this manner 
is very uncertain. Table 6 sunomarizes our values for Vhb and 
Veto- 

To convert Veto to Mv(BTO), we adopt the relationship from 
Gratton et al.'s (1997) analysis of the Hipparcos data, 

MviHB) = 0.17([Fe///] + 1.5) + 0.39. (5) 

In choosing this relationship for the magnitude of the HB, we 
are implicitly choosing absolute ages for the globular clusters of 
about 12 Gyrs. As we discuss in §5.3, this affects the precision 
of our results, but it does not affect our conclusions. 

Chaboyer et al. (1996a) provide the conversion between 
Mv(BTO) and age for V,I data for five metalUcities. The rele- 
vant equations for our discussion are 

tg = 70J-48.5Mv(,BTO) + 93M^^sTO) \Fe/m = -2.5 

(6) 



f9 = 7 1 .4 - 47. SMv^bto) + ^■^M^(bto) {.Fe/H] = -2.0 

(7) 

tg = 86.3-53.5My(Bro) + 7.3M2(£TO)- [Fe///] =-1.5 

(8) 

For all these clusters, we used the equation valid for 
[Fe/Hl=-2.0 to derive ages (Table 6). This is too metal-poor for 
NGC 1466, NGC 2257 and M 3 and too metal-rich for M 92, 
and metallicity does affect the magnitude of the MSTO. To es- 
timate the size of this effect, we used M 3's Mv(bto) to find its 
age using equations (7) and (8), which bracket M 3's [Fe/H]. 
We then assumed a linear relationship between A[Fe/H] and 
A Age and found that we should make M 3's age younger by ~ 
0.6 Gyrs to account for the 0.20 dex difference in metallicity. 
This is stiU consistent with being co-eval with NGC 1466 and 
NGC 2257. A similar calculation using equations (6) and (7) 
suggests that M 92's age should be 0.8 Gyrs older The errors 
shown for A Age in Table 6 include only the observational er- 
rors and not possible errors from incorrect metallicities or dis- 
tances. Hodge 1 1 appears to be 2 Gyrs older by this method, 
but there is a large additional component in the error due to the 
uncertainty in setting the HB level. These calculations confirm 
our result in the previous sections. 

Each method of calculating relative ages suffers from both 
systematic and random uncertainties. However, by using a 
combination of these methods, we can hope to determine the al- 
lowable age differences. Despite potential problems with find- 
ing HB magnitudes, relying on theoretical isochrones for cah- 
bration of age differences in color and magnitude, calculating 
fiducials for clusters with different crowding conditions and us- 
ing the HST calibrations, we consistently find no large age dif- 
ferences between M 92 and M 3 on one hand and NGC 1466, 
NGC 2257 and Hodge 11 on the other. Unless one of these 
methods is found to suffer from a serious flaw, a reasonable 
statement is that the GGC and these LMC clusters have the 
same age ±1.5 Gyr. However, we find that Hodge 11 con- 
sistently looks a little older. This may be due to an actual age 
difference, but the Ay measurement in particular could be due 
to the mixing of helium into the envelope, as discussed above. 
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Table 6 
Veto QUANTITIES 



Cluster 


(V-r,TO 


Veto 


Vhb 


9lhead Mv(eto) 


A Age (Gyrs) 


M3 


0.595 


18.65 ±0.05 


15.75 ±0.05 


3.26 ±0.07 


0.0±0.8 


M92 


0.553 


18.14±0.05 


15.18±0.08 


3.22 ±0.09 


-0.4 ±0.8 


NGC 1466 


0.623 


22.33 ±0.03 


19.32±0.01 


3.34 ±0.03 


0.8 ±0.4 


NGC 2257 


0.610 


22.04 ±0.03 


19.10±0.01 


3.27 ±0.03 


0.1±0.4 


Hodge 11 


0.628 


22.25 ± 0.03 


19.11±0.1 


3.44 ±0.10 


2.0 ±1.5 



The extremely blue HB of Hodge 1 1 , bluer than the more metal- 
poor M 92, would agree with either an age difference or heUum- 
mixing. 

5.5. Distances 

The distance to the LMC is an important and controversial 
quantity (see e.g., Westerlund 1997 for a review). The LMC 
provides a testing ground for consistency between the Cepheid 
and RR Lyrae distance scales. Recent estimates for the LMC 
distance modulus have ranged from (m-M)o = 18.05, based on 
observations of "red clump" stars (Stanek, Zaritsky & Harris, 
1998), to (m-M)o = 18.7, based on Cepheid properties (Feast & 
Catchpole 1997). Femley et al. (1998) used Hipparcos proper 
motions for Galactic RR Lyrae stars and slalislical parallax to 
derive (m-M)o = 18.26. The most recent estimate based on the 
illumination of the SN1987a ring is (m-M)o < 18.44 (Gould 
& Uza 1998). Generally, the RR-Lyrae-based distance to the 
LMC has been smaller than that set by Cepheid observations 
(see e.g.. Walker 1992a). Adding to the confusion is the possi- 
bihty that RR-Lyrae stars in the LMC may have systematically 
different luminosities than Galactic RR Lyraes (Walker 1992a; 
van den Bergh 1995). 

We can use our data to examine these issues by determin- 
ing the relative distance moduli between the LMC clusters and 
M 92 via (1) matching HB levels and (2) matching the un- 
evolved main-sequence (MS) positions. Our data extend far 
enough down the MS that we can use method (2). Because of 
the steep slope of the main sequence, the color offset due to 
reddening and systematic errors in V-I must be removed before 
the MSs can be matched vertically. We do this by shifting the 
RGBs to their expected relative positions and then adjusting the 
distance moduli to match the MSs. This is an iterative process, 
and not as precise as knowing the reddening separately. We find 
that our range of acceptable fits leads to errors in the distance 
moduU of ~ 0.05 mag. As can be seen in Figures 10 through 
12, we would have derived identical relative distances if we 
had matched the HBs of the clusters. This latter observation 
strongly suggests that there is no large difference (~ 0.3 mag 
) in the luminosity of HB stars between the LMC clusters and 
M 92. We then corrected these distance moduli for the differ- 
ential extinction between M 92 and these clusters. We adopted 
the reddenings of Walker (1992a) for our LMC clusters and 
E{B-V) = 0.02mag for M 92 (Harris 1996). These true dis- 
tance moduli between M 92 and the LMC clusters are listed in 
Table 7 as Hm92.gc- 

As discussed earlier, our HST zeropoints have errors ~ 
0.04mag, due to uncertain aperture and CTE corrections and 
the accuracy of the Holtzman et al. calibration. To provide a 
comparison with a ground-based reference, in Table 7 we in- 



clude our HB magnitudes and Walker's < Mv(RR) > corrected 
by 0.04 mag for the offset between RR Lyrae magnitudes and 
HB level (Gratton et al. 1997). Added to this photometric 
uncertainty are our uncertainties in the MS -fitting and in the 
differential reddening. If we take 0.03 mag as an estimate of 
the uncertainty in the differential reddening, then our error per 
cluster is 0.07 mag. 

We prefer the Pont et al. (1998) M 92 distance as it is based 
on a larger, better-selected sample of subdwarfs. The distance 
moduli for M92 and our clusters are listed in Table 7 as Hq^cc- 
Using this, our LMC distance modulus from the average of our 
clusters' distances is m-M~ 18.49. To illustrate the systematic 
uncertainty still remaining because of the dispute over the local 
distance scale, we include in Table 7 the distance moduli we 
derive using the Gratton et al. (1997) distance to M 92. As for 
our measuring error, the MS-fitting error and at least half of the 
zeropoint error are random. The systematic errors, such as the 
errors in the gain ratios of the chips are mitigated because the 
center of the cluster is not always on the same chip. Our clus- 
ters are also not at the center of the LMC. We have attempted to 
correct for this by assuming that the clusters lie in the HI disk 
(inclination 29°, PA of nodes -9°) as argued by Schommer et 
al. (1992). Our derived distances for the LMC center are in- 
cluded in Table 7 ijiQ^iMc)- This changes our LMC distance 
modulus only slightly to 18.52mag. In fact, if the clusters he 
in a disk, then our distance measurement is relatively inmiune 
to changes in the assumed inclination and PA of the nodes be- 
cause NGC 1466 and NGC 2257 have about the same projected 
distance from the LMC center but are on opposite sides of the 
LMC center, while Hodge 1 1 is fairly close to the LMC center. 
For our final error budget, we include the 0.07 mag per cluster 
divided by \/3 plus 0.03 mag error in our averaged number due 
to uncertain geometrical corrections. Also, Pont et al. quote 
0.08 mag as their error in the distance modulus to M92. Our 
final (m-M)o is 18.46 ± 0.09 mag. 

5.6. Blue Stragglers 

All three of our clusters have a prominent blue- straggler star 
(BSS) sequence. To see how the frequency of BSS in these 
clusters compares with that in GGCs, we used the specific fre- 
quency, Fbss, defined by Bolte, Hesser & Stetson (1993) to be 
the ratio of the number of BSS to the number of stars with 
V < Vhb + 2. Before we can count stars in either category, we 
need to make two corrections. First, even the 260s frames were 
long enough to saturate the brightest RGB and AGB stars in the 
clusters and therefore they do not appear on our CMDs. How- 
ever, by looking at our images and star hsts, it is easy to find 
the number of bright stars (~ 15) that are not on the Ust (Table 
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Table 7 
Distance Moduli 



Cluster 


Av 


IJ-M91,GC 




l-'-Q^LMC 




Vhb 










Pont 


Pont 


Gratton 




Walker 


M92 


0.06 




14.61 




14.74 


15.18 




NGC 1466 


0.29 


3.97 


18.58 


18.42 


18.71 


19.32 


19.29 


NGC 2257 


0.13 


3.78 


18.39 


18.55 


18.52 


19.10 


18.99 


Hodge 11 


0.25 


3.89 


18.50 


18.59 


18.63 


19.11 





Table 8 

Stellar Population Statistics 



Cluster 


# "BS" 


# "BS" 


# RGB/HB/AGB 


# Saturated 


# RGB/HB/AGB 


Ratio 




"cluster" 


"field" 


{V<Vhb + 2) 


colheadStars 


{V<Vhb + 2) 


ClusterArea 








"cluster" 




"fie Id" 


FieldArea 


NGC 1466 


73 





721 


17 


26 


1.41 


NGC 2257 


67 


5 


422 


16 


44 


1.64 



8). 

The larger correction is due to field-star contamination in the 
BSS region. The analysis of Hodge 1 1 is severely hampered 
by this, and we will not determine F^ss for this cluster. For the 
"BSS region" in the CMD, we adopted a box with magnitude 
boundaries at V+o.os -1-0 and V+o.os -3.0 and color boundaries 
at V-I MSTO-^-^ and V-I msto^^-(> (see Figures 15-16). These 
boundaries are somewhat arbitrary and are on the conservative 
side (i.e., they minimize contamination from the MSTO, but 
also may miss some BSS). For NGC 1466 and NGC 2257, our 
total sample consists of the stars from all the chips, with the 
shifts discussed in §4 applied (this step makes very little differ- 
ence). To estimate the field contamination, we divided our sam- 
ple at r = 80". The outer stars are not a true field population, 
but will provide an upper limit on the number of interlopers we 
might expect. 

In Figures 15 and 16, we show the upper CMD of both 
the "field" and the "cluster" populations with the BSS region 
marked. The number of stars in each region is included in Ta- 
ble 8. (We note that all of the saturated stars belong to the 
"cluster" sample.) NGC 2257 and NGC 1466 have Fbss of 
0.14 and 0.10, respectively. For these numbers, we have only 
subtracted our "field" contamination from the blue straggler re- 
gion. We have not considered contamination from field stars in 
the RGB-hAGB+HB region. Examination of the CMDs of stars 
further away from the cluster in Walker (1992b) and Testa et al. 
(1995) show that the number of field stars is about equal in both 
regions. Since we have so many bright RGBh-AGBh-HB stars 
compared to blue stragglers, possible contamination in this re- 
gion will not affect our results significantly. In Table 9 we Ust 
the Fbss of some GGCs. In this group we include new numbers 
for Eridanus, Pal 3, Pal 4, and NGC 2419 based on the HST 
observations of Stetson et al. (1999) and Harris et al. (1997) 
and using our definition of BSSs. 

6. DISCUSSION 



6.1. Blue Stragglers: Fbss 

The globular clusters listed in Table 9 have remarkably sim- 
ilar Fbss ■ There are, of course, systematic differences among 
the studies in their completeness, in crowding effects and in 
the blue straggler definition. However, while there may be fac- 
tors of up to two obscured by such problems, factors of ten are 
very unlikely, based on comparing studies of the same clus- 
ter with different completeness limits and crowding conditions. 
We have chosen mainly other HST studies to include in Table 
9, which also makes the sample more homogeneous. The rela- 
tively narrow range in Fbss, compared to the 10"* range in cluster 
central densities, was noted by Sosin & King (1995). They ar- 
gued that as the cluster density changes, different mechanisms 
{e.g., stellar collisions, tidal capture, merging of primordial bi- 
naries) become important for the formation of blue stragglers. 
Although the efficiency of the various mechanisms varies from 
cluster to cluster, the total production rate stays approximately 
constant. Our new data strengthen this point. Eridanus has a 
Fbss at the high end of the range and close to that of M 30. 
M 30 is a post-core-collapse cluster with a central mass density 
log(po)= 5.26. Eridanus, on the other hand, is a sparse cluster 
at large Galactic radius. Its \og{po) is of order 0.5 (Webbink 
1985). The blue straggler specific frequencies for NGC 1466 
and NGC 2257 show that clusters belonging to other galaxies 
have the same constraints on the rate of blue straggler formation 
as those of the Galaxy. 

6.2. Ages and Age Distribution 

6.2. 1 . Previous Results 

The principle result of this study is that these three LMC 
clusters are the same age as M 92 and M 3 to a precision of 1 .5 
Gyr. Mighell et al. (1996) previously reported that Hodge 11 
and M 92 have the same age to within 10 to 21%. Testa et 
al. (1995), based on a comparison of fiducial sequences and 
Ay values, found that NGC 2257 was about 2 Gyrs younger 
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Fig. 15. — The upper part of the CMD of NGC 1466 with the blue straggler box marked. The upper panel shows the stars with r < 80" and the bottom with r > 
80". There is very little field contamination evident in the "field" plot. 



Table 9 
Blue Straggler Frequencies 



Cluster 


Fbss 


Survey Region 


C 


Source for Fbss 






Milky Way Globulars 




M30 


0.19 ±0.04 


r < 10 r core 


2.40 


YannyefaZ. 1994 


NGC 6624 


0.13 ±0.14 


r<l .Orcore 


2.15 


Sosin & King 1995 


47Tuc 


0.07 ± 0.01 


r < 3. Orcore 


2.08 


Guhathakurta a/. 1992 


M3 


0.09 ± 0.02 


r < O.lrcore 


1.89 


Guhathakurta ef flZ. 1994 


M 13 


0.04-0.09 


r ^ 1 fcore 


1.44 


Cohen ef a/. 1997 


NGC 2419 


0.06 ±0.005 


r ^ 8 .2r core 


1.40 


Harris ef a/. 1997 


Eridanus 


0.18 ±0.05 


r< SArcore 


1.14 


Stetson era/. 1999 


Pal 3 


0.15 ± 0.05 


' _ ^•■^1 core 


0.96 


Stetson a/. 1999 


Pal 4 


0.12 ±0.04 


^ ^ 4.9r core 


0.78 


Stetson efflZ. 1999 


NGC 5053 


0.14 ±0.03 


r < 1 .Orcore 


0.77 


Nemec & Cohen 1989 






LMC Globulars 




NGC 1466 


0.10 ±0.01 


r < SArcore 


1.42 


this paper 


NGC 2257 


0.14 ±0.02 


r <3 2r 

' _^ core 


1.10 


this paper 



Note. — Structural parameters for GGCs from Trager et al (1993) and for LMC 
from Mateo (1987). 
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Fig. 16. — The upper part of the CMD of NGC 2257 with the blue straggler box marked. The upper panel shows the stars with r < 80" and the bottom with r > 
80". While there are some stars in the blue straggler box in the outer sample, the overall diagram shows very Uttle evidence for field contamination. Therefore our 
derived Fbss is most likely a lower Umit. 



18 



OLD LMC CLUSTERS 



than NGC 5897 and M 3. In §5.4, we found that a modified 
AV method gives similar ages for M 3 and NGC 2257. Part of 
the discrepancy lies in the smaller AV value of 3.45 we find for 
M 3 (Johnson & Bolte 1998). We note that Stetson et al. (1998) 
found a AV of 3.46 for M3. Their sample included the Johnson 
& Bolte stars, but was about twice as large, including having 
more HE stars. 

There are two other major studies of the ages of other old 
clusters in the LMC. Brocato et al. (1996) argued on the ba- 
sis of AV measurements that NGC 1786, NGC 2210 and NGC 
1841 have the same age as the mean age of the GGCs, although 
their precision was about ± 3 Gyrs. Olsen et al. (1998) (098) 
found that the ages of four old LMC clusters were the same 
as M 5 or M 55. One cluster, NGC 1835, was 2 Gyrs older 
than M 3. Depending on the age difference between M 3 and 
M 5, this could mean than NGC 1835 was coeval with the other 
098 clusters. Another significant result from their paper was 
the determination of the metallicities of their clusters from the 
morphology of the RGB. They found metallicities about 0.2 dex 
higher than those found by Olszewski et al. (1991). 

Before we can put all the above results together, including 
ours, we need to consider possible age differences between the 
GGC samples used to compare with the LMC ones. Lee et al. 
(1994) argued on the basis of HE morphology that GGCs with 
Rqc < 8 kpc are on average 2 Gyrs older than those between 8 
and 40 kpc. M 5 (6.4 kpc) and M 55 (4.2 kpc) fall into the inner 
group while M 3 (8.5 kpc) and M 92 (9.1 kpc) are members of 
the outer, possibly younger halo (although note that the orbit of 
M 5 carries it into the outer halo). Without resolving the issue 
of the chronology of GGC formation, we will discuss results for 
the GGCs listed above. Table 10 summarizes the recent relative 
ages for these four clusters obtained from the literature. 

It is apparent that widely varying opinions exist on the age 
spread in the Galactic halo. Part of this is due to observational 
uncertainties. Our value for the UB magnitude of M 92 (15.18) 
is considerably fainter than the 14.96 used by Chaboyer et al. 
(1996b) and Gratton et al. (1997). If our value is correct, then 
the age gap between M 92 and the other clusters will be re- 
duced. However, although uncertainties still exist, the outer 
halo clusters M 92 and M 3 are apparently not younger than the 
inner halo M 5 and M 55. We can directly compare our LMC 
fiducials to the inner-halo GGC M 55, which has a metallicity, 
|Fe/H|= -1.85 (ZW84), about the same as our LMC clusters. 
G. Mandushev kindly sent us the V, V-I photometry of Man- 
dushev et al. (1996). We choose the NGC 2257 fiducial as 
representative of our LMC clusters. Figure 17 shows that the 
two cluster fiducials match very closely throughout the CMD. 
As is the case for the LMC-M 92 comparison, there is at most a 
small, < 1 .5 Gyr, age difference between NGC 2257 and M 55. 

6.2.2. Implications for the Formation of the Local Group 

The LMC has been suggested as the prototype of a source 
for present-day GGC. One of the scenarios advocated by Searle 
& Zinn (1978) was that clusters formed in "fragments", per- 
haps small gas-rich irregulars, that evolved and interacted over 
the span of several billion years. In the particular case of the 
LMC, Zinn (1980b, 1993) has suggested that satellites Uke the 
LMC contributed their GCs to the MW halo. Previous observa- 
tions of the LMC clusters revealed that they are similar to the 
GGC in many respects (see Suntzeff et al. 1992). Suntzeff et 
al. pointed out that the GGC system could have been made by 
merging ~ 8 LMC-like objects. Van den Bergh (1996) found, 
on the other hand, that the metallicities of the old LMC clusters 



in the Suntzeff et al. sample were, on average, lower than those 
for clusters in the MW's outer halo. The increased metallicity 
estimates by 098 for five of the old LMC clusters would make 
this discrepancy smaller. Van den Bergh & Morbey (1984) and 
others also argued that since the ellipticity of the LMC clus- 
ters is greater than in the MW clusters, it was unlikely that the 
MW had captured LMC clusters. The LMC clusters have a 
larger average elUpticity regardless of the age range considered, 
although at the time, accurate ages were not available. Here 
we consider only those LMC clusters that have been confirmed 
as old and that have ellipticity measurements by Frenk & Fall 
(1982) and those MW clusters that are not heavily reddened 
(E(fi-y)< 0.32) (Harris 1992). A two-sided K-S test can re- 
ject only at the 77% level the hyphothesis that these two groups 
were drawn from the same distribution. Goodwin (1997) notes 
that cluster elhpticity can depend more strongly on the tidal 
field of the parent galaxy than it does on intrinsic cluster prop- 
erties. For the specific case of the three LMC clusters consid- 
ered here, each of them has a small eccentricity, consistent with 
the distribution seen in the GGCs. 

Another similarity between the old clusters of the Galaxy and 
LMC is in the distribution of HB morphology with [Fe/H]. In 
the outer halo of the Galaxy, the HB morphology of GGCs is 
on average too red for the cluster metallicity (the "young halo" 
of Lee et al. 1994). Zinn (1993) showed that the best match 
to the outer halo GGC properties in terms of [Fe/H] and HB 
morphology were the clusters of the LMC and SMC. (In fact, 
the SMC globular NGC 121 was the first "second parameter" 
globular cluster noticed (van den Bergh 1967)). He suggested 
that the majority of the inner halo clusters formed in the over- 
all collapse of the Galaxy, but that the outer halo was mostly 
accreted from relatively large satellites. 

The new HST results allow us to make a few comments about 
the scenario in which LMC-Uke galaxies are absorbed by the 
Galaxy and contribute their stars to the Galactic halo. On the 
basis of age, the old LMC clusters would be indistinguishable 
from existing GGCs, at least our comparison clusters in the in- 
termediate halo. With the recent age estimates and, for the 098 
clusters, revisions in [Fe/H] estimates, it is also possible to im- 
prove on the [Fe/H] - HB morphology comparison between the 
LMC clusters and different GGC populations. In Figure 18, 
we plot our version of HB-type versus [Fe/H] originally shown 
in Lee (1990). The 098 clusters at their revised metallicities 
show nice agreement with the iimer halo clusters of the Galaxy, 
while our clusters fall among the outer halo clusters, as orig- 
inally noted by Zinn (1993). Since the 098 clusters were not 
well-studied before HST and had no HB-type, they were not 
included in the analysis of Zinn (1993). However, with the new 
data, we can see that if the LMC contributes all of its GCs to 
the MW halo, it will contribute clusters like those in the inner 
halo as well, in about equal proportions. 

There is another interesting, though tentative point. If we 
accept the 098 metallicities, then, as is the case for the GGC, 
the LMC clusters also separate into two groups in the HB-type 
versus [Fe/H] plot. One group, composed of the 098 clusters, 
has higher [Fe/H]'s but not redder HBs than the other, which 
includes all the clusters studied by Walker (1992b), including 
NGC 1466, NGC 2257 and Hodge 11. The LMC groups also 
mimic the GGC in their spatial distribution differences. The 
098 clusters are concentrated toward the center of the LMC. 
They all have projected angular separations < 2.5° from the 
LMC center, whereas the Walker sample has no clusters closer 
than 2.5° from the center. Based on our analysis in §6.2.1, 
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Table 10 

A Ages (Gyrs) for Comparison GGCs 



Source 


M92 


M55 


M3 


M5 


Method 


Richer era/. (1996) 


0.0 


0.5 


-1.7 




5iV-D 


Chaboyer ef aZ. (1996b) 


0.0 


-2.6 


-4.9 


-4.3 


AV 


Grattoneffl/. (1997) 


0.0 






-2.9 


MS fitting 


Buonanno el al. (1998) 


0.0 


1.1 


-0.2 


-0.2 


Ay. MV-I) 



14 



16 



> 

IB 



20 



0.5 1 1.5 

V-I 

Fig. 17. — NGC 2257's fiducial compared to the M 55 data of Mandushev et al. (1996). The overall agreement between these two clusters of comparable 
metalUcity shows that there is not a large age difference between them. 
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Fig. 18. — HB-type versus [Fe/H] for LMC and Galactic globular clusters. The HB-type is defined to be the number of blue stars minus the number of red stars 
divided by the total number of stars on the HB, including variables. The Galactic clusters have been divided into two groups, as suggested by Lee et al. (1994). 
Those MW clusters with R < 8 kpc tend to have bluer HBs than similar metallicity ones with R > 8 kpc. The LMC clusters display the same trend when the inner 
and outer clusters are compared. The data for inner LMC clusters are from Olsen et al. (1998), while the data for the outer ones are from Walker (1992b). The data 
for the MW come from the compUlation in Lee et al. (1994). We have marked with a dotted circle those MW clusters that are most Ukely truly young clusters, based 
on their tumoff magnitudes. They are concentrated at even redder HBs than the LMC clusters. 



where we discussed the age range among all the LMC clus- 
ters studied with HST, it is possible that age in this case in not 
the second parameter, a point also made by Da Costa (1999). 
A picture where the LMC builds up its outer GCs by accreting 
satellites is also less appeahng, because of the mass of the LMC 
is not much larger than the the mass of the smallest satellites 
known to have their own globular cluster systems. An addi- 
tional clue to the formation of the LMC cluster system comes 
from the possible radial abundance gradient, if the 098 metal- 
Ucities are correct (Da Costa 1999). 

6.2.3. The Big Picture 

Our most fundamental result is that the oldest LMC clus- 
ters are the same age as the old GGC. What does this mean for 
our understanding of galaxy formation? On one extreme, if we 
imagine the Galaxy being made completely via the accumula- 
tion of LMC-sized fragments (and if we accept the argument 
of Stetson et al. (1996) that most GGC are essentially coeval), 
then this coincidence in age would be telling us that the collapse 
epoch of all the original fragments was the same to 1.5 Gyr or 
so. In the opposite extreme, in which only a small fraction of 
the current-day Galaxy was acquired through accumulation of 
small companions, our result is then telUng us that the forma- 
tion epoch for the Milky Way out to the intermediate halo is 
the same to ~ 1.5 Gyr as for the LMC, a galaxy ten times less 
massive. 

We discuss briefly here this result in the context of one (pop- 
ular) galaxy formation scenario - hierarchical clustering (e.g., 
Blumenthal et al. 1984). It is not yet possible to fully simu- 
late the processes that lead to galaxy formation. In particular, 
the question of when globular clusters form within hierarchical 
clustering models is not yet accessible to simulations. Never- 
theless, with the reasonable assumptions that globular clusters 
were formed within larger structures (e.g. a proto-LMC) and 



were among the first objects formed after the gas began to col- 
lapse, we can investigate whether the current hierarchical for- 
mation models can be consistent with our "fossil record" data. 
We note that as the models become more complete and sophisti- 
cated, this discussion may become rapidly dated. Nevertheless, 
as the data for globular cluster ages continue to improve, they 
will become increasingly important in testing models such as 
these. 

As an example of what hierarchical clustering models predict 
for the formation epoch of LMC-like galaxies, we examined the 
properties of galaxies in 100 "Milky Way" halos calculated by 
the semi- analytic method of Somerville & Kolatt (1999) and 
Someville & Primack (1999) . These semi-analytic models are 
a powerful way to study galaxy formation by combining merg- 
ing histories for dark-matter halos based on extended Press- 
Schechter theory (Bower 1991; Bond et al. 1991) with simple 
prescriptions for such physical processes as gas cooling, star 
formation and supernova feedback. The advantage of using this 
method is that the formation and evolution of many systems in 
many different cosmologies can be calculated in a reasonable 
period of time. 

Placing GC formation in this scenario is difficult, because 
there is no agreed-on model of globular cluster formation and 
these models are not hydrodynamical calculations that can 
probe globular cluster formation. However, since the globular 
clusters are sUghtly metal-enriched and since models of glob- 
ular cluster formation generally start with cold molecular gas 
(e.g., Harris and Pudritz 1994; Elmegreen & Efremov 1997), 
we can place the epoch of formation just after the gas can cool 
and start to form stars. When this happens depends on the mass 
of the dark matter halo. Because of the meta-galactic UV field, 
gas is prevented from cooling until its dark matter halo has 
reached the critical mass when its gas is dense enough to be 
self-shielded from ionizing photons. Simulations have shown 



JOHNSON ET AL. 



21 




02468 10 02468 10 
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Fig. 19. — A histogram of the time lapse between the collapse of the first progenitor of the parent galaxy and of each of its bright (Mg-Slogh < -15) satellites 
for our "Millcy Way" halos in (a) a standard CDM universe and (b) a ACDM universe. At the level that we can determine age differences, these two distributions 
are essentially the same. 



that while this critical mass depends on redshift, it is not sen- 
sitive to such details as the shape of the UV-spectrum and the 
inclusion of star formation, especially at the level of precision 
we require {e.g., Weinberg, Hernquist & Katz 1997). 

For each present-day galaxy in our "Milky Way" halo cata- 
log, we have the redshift (zcoii) when any of the halos in that 
particular galaxy's merging tree reached that critical mass in 
dark matter (the "first progenitor"). We assume that this marks 
the initial stage of globular cluster formation. We confine our 
results here to those satellites with Mb < 51og(h)- 15^, which 
includes the LMC and SMC, but not the dwarf spheroidals. 

To illustrate the range of separation in time of this epoch, we 
calculated the difference in Zcoii between the Milky- Way type 
galaxy and each of its satellites in the 100 model halos. The 
results are plotted in Figures 19a-b for two different cosmolo- 
gies, standard cold dark matter (SCDM) with n,„atter=i-0 and 
H()=50 km s"' Mpc"' and cold dark matter with a cosmologi- 
cal constant (ACDM) with = 0.7, i^matter = 0.3, and Ho=70 
km Mpc~' . While there is a tail to large age differences, 
there is a strong peak at small age differences in both cosmolo- 
gies. We cannot measure age differences much smaller than 1 
Gyr. Therefore, if we assume that the initial epoch of globu- 
lar cluster formation in a galaxy happened at about the same 
time as the collapse of the first progenitor of that galaxy, then 
our result that the oldest LMC and the Galactic clusters formed 
at the same time ±1.5 Gyr is completely consistent with cur- 
rent hierarchical clustering models. As Harris et al. (1997) 
pointed out, this apparent synchronicity between the LMC and 
Galaxy, and within the outer halo of the Galaxy, is due to a 
number of individual dark mailer halos approximately simulta- 
neously reaching a critical mass so that star formation can start. 
The SMC provides a counter-example as it apparenlly began 
its cluster formation ~ 2 Gyr later {e.g., Stryker, Da Costa, & 
Mould 1985; Mighell, Sarajedini, & French 1998). In our hier- 
archial clustering models, it would be in the tail of the "Milky 
Way"-satelhte age distribution. 

We conclude that hierarchial clustering models, at least as 
described by current semi-analytical codes, are consistent with 
both constraints mentioned at the beginning of this section, 
though it is sUghtly harder to satisfy the first constraint for a 



ACDM case. While many LMC-size galaxies form at about the 
same time, there is dispersion in Zmii, making it less likely to 
have all the original fragments collapse at the same time. On 
the other hand, the mass difference between the Milky Way and 
the LMC is not enough for it to be unusual for the collapse of 
the first progenitor to happen at the same time for both galaxies. 
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